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Abstract 

We examine a ^^P donor electron spin in a Si crystal to be used for the purposes of quantum 
computation. The interaction with an uncontrolled system of ^^Si nuclear spins influences the 
electron spin dynamics appreciably. The hyperfine field at the ^^Si nuclei positions is non-collinear 
with the external magnetic field. Quantum operations with the electron wave function, i.e. using 
magnetic field pulses or electrical gates, change the orientation of hyperfine field and disturb the 
nuclear spin system. This disturbance produces a deviation of the electron spin qubit from an ideal 
state, at a short time scale in comparison with the nuclear spin diffusion time. For Hgxt ~ 9 T, the 
estimated error rate is comparable to the threshold value required by the quantum error correction 
algorithms. The rate is lower at higher external magnetic fields. 



*Electronic address: |ssaykin@clarkson.edti| also at Physics Department, Kazan State University, Kazan 
420008, Russia 

'''Electronic address: lleonid© clarkson . edul 



1 



The interest in the evolution of a single electron spin interacting with a system of uncon- 



trolled nuclear spins 
qubit manipulations 




3ath) arose 



rom proposed models of spin qubits and qubit- 
8|. The electrons confined in a quantum dot or 



in a hydrogen state around a shallow donor impurity are influenced by many nuclear spins 



through the 
decoherence 



lyperfine dipole-dipole interaction. At low temperatures, where spin-phonon 
9j is ineffective, the dominant mechanism for the electron spin decoherence is 
via this interaction 13^ . Additional error in a qubit state appears from non-ideal 

quantum gate operations [1J|. In this work we consider the effect of the nuclear spin bath 
on results of the quantum operations. 

The approximation of short bath correlation times and weak system-bath interaction, 
which leads to a Markovian-like equation for relaxation and decoherence, is hardly applicable 
for quantum computation. In that case the evolution of qubit at short times is important jl5|. 
Different terms of system-bath interaction and the bath dynamics influence the evolution of 
a single qubit at different time scales. In particular, for the spin qubit, interacting with a 
spin bath, several time scales can be defined 

The irreversible dynamics of the nuclear spin bath determines the longest time scale. At 
low temperatures this can be explained in terms of spin diffusion \v\. The characteristic 
diffusion time is determined by the fiip-fiop transitions (/|/^ + 1-1%) for two neighbour- 
ing nuclear spins. These processes are much slower in comparison with the electron spin 
dynamics. For shorter time intervals, the electron spin sees a "frozen" system of nuclear 
spins. The diagonal component, Szlz-, of the hyperfine interaction with a "frozen" spin sys- 
tem leads to the constant shift of the electron Larmor frequency. It is possible to eliminate 
this in quantum algorithms by various additional correction operations (in analogy with a 
7r/2 — TT pulse sequence), or by the preliminary measuring of the "nuclear magnetic field" 
at the localized electron. The single transitions with the electron spin flip (S'+Jq, + S'_/q,, 
a = +,—,z) are suppressed by energy conservation, unless the external magnetic field is 
sufficiently low [l^ll^l- However for nuclear spin flips, this suppression is not as strong, due 
to the small Zeeman splitting of nuclear levels. 

The effective magnetic field at each nuclear position is formed by the external and hyper- 
fine fields. It is non-colinear to the external field, due to the hyperfine terms, Szl+ + Szl^- 
Quantum operations with the electron spin reorient this field, and cause the nuclear spin 
bath to oscillate near the nuclear Larmor frequency. These oscillations can be appreciably 
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faster compared to nuclear spin diffusion. As a result, this produces a deviation of the elec- 
tron spin qubit from an ideal state. We study this effect in detail in the case of a ^^P donor 
electron spin qubit system. 

If, initially, at the time t = 0~ the qubit is uncoupled with the nuclear spin bath, and 
the bath is in an equilibrium state, Po-(0^)pbat/i(0~), after very fast quantum operation the 
system-bath state will be Pa{0'^)pbath{0^)- We assume here, that the qubit operation does 
not change the bath state significantly. The qubit density matrix at the time t is 

p,{t) = Tr,,,, {e-i^V.(0+)p6a*.(0-)e^^*} . (1) 

At short times, in comparison with the nuclear spin diffusion, we neglect by the nuclear 
dipole-dipole interaction. The magnetic Hamiltonian of a single electron spin localized at a 
phosphorus donor in a ^^Si^; ^^Sii_^. crystal will be 

H = (7/iBHS + aoSIo - 7p^HIo + Yl (SA(Rfc)Ifc - 7si^HIfc) , (2) 

k 

where Iq is a donor nuclear spin and 1^ are ^^Si nuclear spins. The electron spin interaction 
with the phosphorus nucleus, the second term on the right-hand side of Eq. Q, is isotropic. 
The tensor, A(Rfc), describes the hyperfine interaction of an electron spin with a nuclear 
spin of ^^Si at a lattice position R^. It consists of isotropic and anisotropic parts as 

g 

where we use the free electron (yf-factor Values of hyperfine constants have been mea- 
sured by G. Feher [20] and later by E. B. Hale and R. L. Mieher Q in the ENDOR ex- 
periments for different positions of nuclear spins. At some lattice positions the off-diagonal 
elements of tensor A(Rfc), Eq. (jHl), are comparable with the diagonal elements. Table HI 
Though the effective mass theory desribes symmetry of the hyperfine tensors incorrectly 



the reasonable results have been obtained in the Bloch representations 
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22 1 or in 



Wannier representations 



23| of a '^^P donor electron wave function. 
To calculate an evolution of the electron spin density matrix with the Hamiltonian, Eq. 
0, we use the approach derived by W.B. Mims (2J] for the modeling of electron spin echo 
envelope modulation functions (ESEEM). In a high magnetic field, we assume that the effect 
of the non-diagonal terms (S'+Jq, + S'_/q, a = +,—, z) in Eq. (0) is small. These transitions 
are forbidden in the first order of perturbation and inversely proportional to the electron 



n 



TABLE I: Constants of the hyperfine interaction from the ENDOR experiments |21l | and numbers 

of equivalent positions, n, near the ^^P donor. 

Shell |R| (A) (KHz) A,^ (KHz) 

E (1 1 1) 2.35 540 700 4 

B (4 4 0) 7.68 4474 -39.8 12 

H (4 4 -4) 9.41 1378 -50.6 12 

R (7 7 1) 13.51 779 19.4 12 



Zeeman level splitting in the second order. The truncated Hamiltonian can be written in 
the form 



(4) 



It is diagonal for the electron spin but non-diagonal for spins of ^^Si nuclei. In the interaction 
representation with 

= gfi^RS, + a^SJl - 7pmJ°, (5) 

the Hamiltonian, Eq. (jH), can be diagonalized for different eigenvalues of Sz using the 
rotation operator 

S ^J^^^Vk{^s)I'^^^e^I^.^ (6) 

it 



where tan6'fe = A^j^(Rfc)/v4^^(Rfc), A^(Rfc) = ^ Al^{Rk) + Aly{Rk) and tanv3fc(mj = 
msA^{Yik)/{Azz{^k)ms - 7si^H). 

We consider a 7r/2 magnetic pulse as the single qubit operation. Op, and assume that 
initially {t = 0^) the nuclear system is in equilibrium and unpolarized. The electron spin 
density matrix in the interaction picture, Eq. (0), at time t is, 

(P<xW)m.< = 11 jcOS j COS 1 j 

+ sin^ ^ MO-Mms) )^ I^ MO+Mms) ;)^ | (,^(o+))^^^^ , (7) 

where the effective oscillation frequency of the nuclear spin located at the lattice site k with 
the electron spin projection is 



Uk{ms) = l/h^J{Azz{Rk)ms - 7si^H)2 + (A^(R^)^j2_ 
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In a high magnetic field, where 

H> |A,,(R,,)/(27si/i)|, (9) 
this effect can be treated as a perturbation, and Eq. ((Tj) is simphfied to 

The reorientation of the effective magnetic field at each nuclear spin position changes the 
static component of the hyperfine field at the bounded electron, the first term in Eq. (fTUIl . 
and admixes oscillations near the ^^Si nuclear Larmor frequency, the second term in Eq. pup. 
The amplitude and width of the admixed band are proportional to H^^. This indicates that 
the effect can be reduced to any expected value by tuning of the external magnetic field. 
In the case of a polarized nuclear spin bath, the polarization pk for each nuclear spin, in 
analogy with a spin temperature, should be used. Therefore, Eq. ((Tj) transforms to 
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In a particular qubit implementation, the qubit density matrix will depend on the position 
and polarization of each ^^Si spin. We estimate the value of this field reorientation effect 
for a ^^P electron spin interacting with a single ^^Si nuclear spin placed at one of the lattice 
sites from Table H] in an external magnetic field, H = 1 T. The magnetic isotope ^^Si, which 
is located in one of the four, nearest to the donor, lattice positions, experiences the largest 
reorientation of the effective magnetic field due to the electron spin manipulation. Fig. ^ If 
the nuclear spin is initially oriented parallel to Hes, Fig. ^a), then after the electron spin 
fiips, it will rotate about new axis of quantization. Fig. ^b). To quantify the infiuence of 
this rotation on the quantum bit evolution, we calculated the norm of deviation of qubit 



state from ideal one 



25| 



\Ap^m = \Mt) - PTim, (12) 



where the deviation of the norm for single qubit is H^H = \/ A11A22 + A12A21. As the ideal 
state we used the spin polarized state, Eq. (fTT|) . without high frequency components. 

The simulated time evolution of the density matrix deviation, after Op = tt/2, is shown 
in Fig. 121 for an initial nuclear spin polarization, pk = 1. For time scales on the order of 100 
ns, the amplitude of the density matrix deviation varies from to 7 x 10"'^, depending on 
the nuclear spin position. 

To perform practically useful large scale quantum calculations, quantum error correction 
procedures should be used. The best known quantum error- correct ion schemes require the 



error rate to be less than one error per sequence of 10^ — 10^ quantum gates 26|, |27|, |28| . Our 
results indicate, for solid state nuclear spin quantum computers, the considered mechanism 
can be crucial, especially for qubits in which the ^^Si isotopes are located at the nearest 
positions to the ^^P. 

For natural silicon crystal, n(^^Si)~ 4.76%, the probability to have a ^^Si in one of the 
four neighbors to the donor position is approximately 18%. Isotopic purification of a silicon 
crystal will reduce the probability to have a "wrong" qubit only. To suppress undesirable 
qubit dynamics, one should use high magnetic fields. It is possible to estimate the value of 
threshold magnetic field, using the || Ap 

maxll ~ H dependence, from -Q — 1 relation, 
where Q ~ 10^ is quality factor required for application of quantum error correction schemes. 
If we assume the effect of a whole nuclear magnetic bath is of the order of the maximal 
single nuclear spin effect, then the threshold value of the external magnetic field should be 
Hth ~ 9 T. 

In the opposite case, to analyze this small effect, one should define experimental conditions 
in which it has considerable influence. To investigate the electron spin echo modulations, 
due to a nuclear spin at the E position (TableHJ, the magnetic field should be, at least, of the 
order of 0.1 T (the electron Zeeman frequency uz ~ 2.8 GHz). Although the ESEEM, 
due to interaction with ^^Si nuclear spins, has been observed on the irradiation- induced 
defects in quartz at the electron Zeeman frequency uz = 1.7 GHz [2^, as far as we know, 
there was no published works on spin echo modulations for the phosphorus donor electron 
spin. 

In the case of GaAs or InGaAs quantum dots the concentration of nuclear magnetic 
isotopes is 100%. The electron spin interacts with a much larger set of a nuclear spins, and 
we expect more complicated dynamics of system. While there are no experimental values 
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of hyperfine constants for GaAs and InGaAs quantum dots, it may be possible to estimate 
them theoretically. We suppose that, in analogy with silicon, calculations in the effective 
mass approximation would be incorrect and one should use, at least, the LCAO model. 

In conclusion, we considered the effect of the non-diagonal terms in the dipole-dipole 
interaction of a ^^P electron spin with a system of ^^Si nuclei in a high magnetic field. 
Electron spin manipulations generate fast oscillations of the nuclear spin bath. Due to this 
effect the electron spin density matrix deviates from an ideal state. In an external magnetic 
field of H = 1 T, the value of deviation is on the order of 7 x 10~^. This is crucial for 
the implementation of practically useful quantum algorithms with electron spin qubit. The 
estimated threshold value of the magnetic field for a fault tolerant ^^P electron spin qubit 
is Hth fti 9 T. 
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FIG. 1: The reorientation of the effective magnetic field at the lattice position (1,1,1) due to the 
electron spin flip. and Hj^ are z and in-plane components of the effective magnetic field Heff at 
the ^^Si position. 
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FIG. 2: The norm of the density matrix deviation from the ideal state as a function of time for 
^^Si isotopes located in different lattice positions, see Tabled (a) position E, (b) position B, (c) 
postion H, (d) position R. The external magnetic field H= 1 T. 
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